Distinguishing separable and entangled states 
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We show how to design famílies of operational criteria that distinguish entangled from separable 
quantum states. The simplest of these tests corresponds to the well-known Peres-Horodecki positive 
partial transpose (PPT) criterion, and the more complicated tests are strictly stronger. The new 
criteria are tractable due to powerful computational and theoretical methods for the class of convex 
optimization problems known as semidefmite programs. We successfully applied the results to many 
low-dimensional states from the literature where the PPT test fails. As a byproduct of the criteria, 
we provide an explicit construction of the corresponding entanglement witnesses. 



Entanglement is one of the most striking features of 
quantum mechanics. Not only is it at the heart of the 
violation of Bell inequalities [1], but it has lately been 
recognized as a very useful resource in the field of quan- 
tum information. Entanglement can be used to perform 
several important tasks such as teleportation, quantum 
key distribution and quantum computation [2], Despite 
its widespread importance, there is not a procedure that 
can teli us whether a given state is entangled or not, and 
considerable effort has been dedicated to this problem in 
recent years [3, 4]. In this letter we apply powerful tools 
of optimization theory for problems known as semidefi- 
nite programs to construct a hierarchy of tests that can 
detect entangled states. 

A bipartite mixed state p is said to be separable [5] (not 
entangled) if it can be written as a convex combination 
of pure product states 

p = Y^Pi\Í>i){Í ) i\®\ ( Í>i)( ( t>i\i (!) 

where \tpi) and \<j>ï) are state- vectors on the spaces Ha 
and Hb of subsystems A and B respectively, and pi > 
0, YliPi = 1- a state admits such a decomposition, 
then it can be created by local operations and classical 
communication by the two parties, and hence it cannot 
be an entangled state. 

Several operational criteria have been proposed to 
identify entangled states. Typically these are based on 
simple properties obeyed by all separable states and are 
thus necessary but not sufficient conditions for separabil- 
ity (although some sufficient conditions for separability 
are known [6] ) . The most famous of these criteria is based 
on the partial transposition and was first introduced by 
Peres [7]. It was shown by the Horodeckis [8] to be both 
necessary and sufficient for separability in H2 <8> H2 and 
H2®Hz- If p has matrix elements Pik,ji = {i\®(k\p\j)Ç$>\l) 
then the partial transpose p TA is defined by p[ k A ^ = Pjk,a ■ 
If a state is separable, then it must have a positive partial 
transpose (PPT). To see this consider the decomposition 
(1) for p. Partial transposition takes \ipi}(ipi\ to \tp*}{tp*\, 
so the result of this operation is another vàlid dcnsity ma- 



trix and must be positive. Thus any state for which p A 
is not positive semidefmite is necessarily entangled. This 
criterion has the advantage of being very easy to check, 
but there are PPT states that are nonetheless entangled 
as was first demonstrated in [9]. 

Our separability criteria will also be based on simple 
computationally checkable properties of separable states. 
Consider the state p defined on Ha®Hb®Ha, given by 

p = ^Pilipi)^ <8> \4>i){<t>i\ ® IV'iX^I· (2) 

Firstly p is an extension of p (that is the partial trace over 
the third party C is equal to p, Trc [p] = p). Secondly 
the state is symmetric under interchanging the two copies 
of "Ha- To put this more formally we define the swap 
operator P such that P\i) <g \k) ® \j) = \j) ® |fc) <g> \i). 
We have P 2 — I, and tt = (7 + P)/2 is a projector onto 
the symmetric subspace. Since Trpn = p the extension p 
only has support on this subspace. Finally the extension 
p is a tripartite separable state. This means that it will 
have positive partial transposes with respect to any of the 
parties, and in particular we have p TA > and p TB > 0. 

We may now formulate an explicit separability crite- 
rion based on the existence of the extension discussed 
above. If the state p on Ha®Hb is separable then there 
is an extension p on Ha ® Hb ® Ha such that Trpn = p, 
p TA > and p TB > 0. Note that the symmetry of the ex- 
tension means that if p TA > then p Tc > 0, so including 
this would not make a stronger test. We may generalize 
this criterion to an arbitrary number of copies of both 
Ha and Hb- If the state p on Ha ® Hb is separable 
then there is an extension p with support only on the 
symmetric subspace of 7í® fc <£>H^ 1 such that p has a pos- 
itive partial transpose for all partitions of the k+l parties 
into two groups. Since the extensions are required to be 
symmetric, it is only necessary to test the possible parti- 
tions into two groups that are not related by permuting 
copies of Ha and Hb- Including testing for positivity of 
the extension itself, there are \(k + 1) (l + 1) /2] distinct 
positivity checks to be satisfied by p. 

These results generate a hierarchy of necessary condi- 
tions for separability. The first is the usual PPT test for 
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a bipartite density matrix p. If the test fails, the state 
is entangled; if the test is passed, the state could be sep- 
arable or entangled. In the latter case we look for an 
extension p of p to three parties such that irpm = p that 
satisfies the PPT test for all possible partial transposes. 
If no such extension exists, then p must be entangled. If 
such an extension is possible, the state could be separa- 
ble or entangled, and we need to consider an extension 
to four parties and so on. 

Each test in this sequence is at least as powerful as the 
previous one. We can see this by showing that if there is a 
PPT extension p n to n parties, then there must be a PPT 
extension p n -i to n — 1 parties. Let p n -\ = Tvx[p n ], 
where X represents one of the copies of A or B. It is 
easy to check that p n -i will inherit from p n the property 
of having its support on the symmetric subspace. Let's 
assume that it is not PPT. Then there is a subset X of the 
parties such that p^Li has a negative eigenvalue, where 
Tx represents the partial transpose with respect to all 
the parties in subset X. Let |e) be the corrcsponding 
eigenvector and let {\i)} be a basis of the system X ovcr 
which the partial trace was performed. Since p n is PPT, 
then (e\(i\pZ x \e)\i) > 0, for all i. Then 

^(e\(i\p^\e)\i) = (e\Tr x [p^]\e)>o. ( 3 ) 

i 

Since X £ X , we can commute the trace and the 
partial transpose, and using p n -i = Trx[/5 n ], we have 
i e \Pn-i\ e ) — 0' which contradiets the fact that e) is an 
eigenvector of p^_\ with negative eigenvalue. 

The problcm of scarching for the extension can be 
solved efficiently, since it can be stated as a particu- 
lar case of the class of convex optimizations known as 
semidefinite programs (SDP) [10]. A SDP corresponds 
to the optimization of a linear function, subject to a lin- 
ear matrix inequality (LMI). A typical SDP will be 

minimize c T x 

subject to F(x) > 0, (4) 

where c is a given vector, x = (xi, . . . , x m ), and F(x) = 
Fq + J2i x iFíi f° r some fixed n-by-n hermitian matrices 
Fj. The inequality in the second line of (4) means that 
the matrix F(x) is positive semidefinite. The vector x 
is the variable over which the minimization is performed. 
In the particular instance in which c = 0, there is no func- 
tion to minimize and the problem reduces to whether or 
not it is possible to find x such that F(x) is positive 
semidefinite. This is termed a feasibility problem. The 
convexity of SDPs has made it possible to develop sophis- 
ticated and reliable analytical and numerical methods for 
them [10]. 

The separability criteria we introduced above may all 
be formulated as semidefinite programs. For brevity we 
will cxplicitly consider only the problem of searching for 
an extension of p to three parties. We will also relax the 



symmctry requirements on the extension p, and we will 
ask only PpP = p. This increases the size of the SDP, but 
simplifies the setup. Let {<jf}i=i t ... i d 2 A A <T f}j=i,...,d 2 B be 
bases for the space of Hermitian matrices that operate 
on Ha and Hb respectively, such that they satisfy 

Tr(o-f af ) = aóij and Ti(af) = 6 il , (5) 

where X stands for A or B, and a is some constant — the 
generators of SU(n) could be used to form such a basis. 
Then we can expand p in the basis {af <S> af}, and write 

P = Eij Pii°i ® of, with Pij = a~ 2 Tr[paf © af}. Wc 
can write the extension p in a similar way 

P = J2pk3i{<rf®<r? ® ef + <?k ® <?? ®°f} + 

ij 

i<k 

+ J2P^ ka k^ a f^ <7 t (6) 

kj 

where we have explicitly used the symmetry between the 
first and third party. We also need to satisfy Trc(f>) = 
p. Using (5), and the fact that the af (S> af form a 
basis of the space of hermitian matrices on Ha'3'Hb, we 
get piji — pij. The remaining components of p will be 
the variables in our SDP. The LMIs come from requiring 
that the state p and its partial transposes be positive 
semidefinite. For example, the condition p > will take 
the form F(x) = F + ^\ x l F i > if we definc 

F = PiJ a i ® of <S> af + 

3 

i=2,j=l 

= of ®af ®af i>2, 
F i3 k = [vf ® af ®af + af ® af ® af) k>i> 2. 

The coeíficients pijk(k ^ l,k > i) play the role of the 
variable x. There are m = {d A A d 2 B — d 2 A d 2 B )/2 compo- 
nents of x, where dj is the dimension of Hi . Each F is a 
square matrix of dimension n = d 2 A ds. Positivity of the 
partial transposes Ta and Tb leads to two more LMIs, 
p TA > and p TB > 0. The F matrices for these two 
LMIs are related to the matrices F^t- by the appropriate 
partial transposition. We can write these three LMIs as 
one, if wc define the matrix G = p © p TA © p TB , so for 
example G = F © F^ A © F^ B (a block-diagonal matrix 
C = A © B is positive semidefinite iff both A and B are 
positive semidefinite) . So the feasibility problem reduces 
to attempting to find pijk(k ^ 1, k > i) with G > 0. In 
fact, the SDP corrcsponding to minimizing t subject to 
Haba + G > is always feasible and performs better 
numerically. A positive optimum gives a value of p* such 
that (l—p)p+plAB l 'dAds is entangled for all < p < p* . 
Looking for an extension on 7Y® fe ® H% 1 is a semidefinite 
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program with m = ( dA+ k k 1 ) 2 { dB+ l 1 ) 2 — d 2 A d 2 B vari- 
ables and a matrix G with \(k + l)(/ + l)/2] blocks of 

dimcnsion at most (^E^" 1 ) 8 ( ) 

Numerical SDP solvers are described in detail in [10]. 
Typically they involve the solution of a series of least 
squares problems each requiring a number of operations 
scaling with problem size as 0(m 2 n 2 ). For SDPs with a 
block structure these break into independent parts each 
with a value of n determined by the block size. The num- 
ber of iterations required is known to scale no worse than 
0(n 1//2 ). Thus for any fixed value of (k, l) the computa- 
tion involved in checking our criteria scales no worse than 
0(d A 3k ^ 2 d B 31 ^ 2 ) which is polynomial in the system size. 

Using the SDP solver SeDuMi [11], wc applicd the first 
criterion (k = 2,1 = 1) to several examples of PPT en- 
tangled states with Úa = 2, ds = 4 or dA = 3,ds = 3. 
On a 500 MHz desktop computer a single state could 
be tested in under a second for dA = 2, ds = 4 and in 
around eight seconds for dA = 3, ds = 3. For the onc 
and two parameter famílies of PPT entangled states de- 
scribed in [3, 9, 12] we performed a systematic search 
of the parameter space, in each case testing hundreds 
or thousands of different states. We checked 4000 ran- 
domly chosen examples of the seven parameter family of 
PPT entangled states states in [13]. We also checked 
the PPT entangled states constructed from unextendiblc 
product bases in [14]. We did not find any PPT entan- 
gled state with an extcnsion of the required form, thus 
vcrifying the entanglcment of all these states. Very close 
to the separable states the test was inconclusive due to 
numerical uncertainties. Uncertainties and one examplc 
are discussed more fully below. 

A very useful property of a SDP, is the existence of 
the dual problem. If a problem can be stated as a SDP 
like (4), usually called the primal problem, then the dual 
problem corresponds to another SDP, that can be written 

maximize —Tt[FqZ] 
subject to Z > 

Tr[F t Z] = a, (7) 

where the matrix Z is hermitian and is the variable over 
which the maximization is performed. For any feasible 
solutions of the primal and dual problems we have 

c T x + Tr[F Z] = Tr[F(x)Z] > 0, (8) 

where the last inequality follows from the fact that both 
F(pí) and Z are positive semidefinite. Then, for the par- 
ticular case of a feasibility problem (c = 0), cquation (8) 
will read Tr[FoZ] > 0. This result can be used to give a 
certificate of infeasibility for the primal problem: if there 
exists Z such that Z > 0, Tr[FiZ] = 0, that satisfies 
Tr[F Z] < 0, then the primal problem must be infeasible. 

In the context of entanglement, the role of the "cer- 
tificate" is played by observables known as entanglement 



witnesses (EW) [8, 15]. An EW for a state p satisfies 

Tr[p sep W] > and Tr[pW] < 0, (9) 

where p sep is any separable state. If our primal SDP 
is infeasible (which means that the state p must be en- 
tangled) , the dual problem provides a certificate of that 
infeasibility that can be used to construct an EW for p. 

First, we note that due to the block diagonal structure 
of the LMI, we can restrict any feasible dual solution Z 
to have the same structure, i.e., Z = Z © Z\ A © Zj B 
where the Zi are operators on TÍa ®Hb® Ha- Then we 
have that Tv[G Z] = Tt[F {Z + Z x + Z 2 )\. We defined 
Fq as a linear function of p so that F = A(p) where A 
is a linear map from Ha <8> Wb to Ha ® Hb ® Ha- Wc 
can now define an operator Z on Ha ®Hb through the 
adjoint map A* such that Z = A*(Z + Z\ + Z 2 ) and 

Tr[pZ] = Tr[A(p)(Z + Z 1 + Z 2 )] = Tt[G Z]. (10) 

If p sep is any separable state, we know that the primal 
problem is feasible (the extension p exists). Then, using 
Tv[G Z] > and (10), wc have Tr[p sep Z] > for any Z 
obtained from a dual feasible solution. For this particular 
problem, if the primal is not feasible (which means p is 
an entangled state), a feasible dual solution Zew that 
satisfies Tt[GoZew] < always exists. Using (10) we 
can see that the corresponding operator Zew satisfies 
Tt[pZew] < which together with Tr{p sep ZEw] > 
means that Zew is an entanglement witness for p. 

In numerical work, if the SDP solver cannot find an 
extcnsion p it construets the matrices Zi. Evaluating 
Tt[pZew] and vcrifying the three positivity conditions 
provides an independent check of the result. Unless this 
check is not conclusive — for example, if Tt[pZew] is not 
significantly different from zero — we are able to defini- 
tively conclude that no p exists. 

If W is an EW, then for any product state \xy) we have 

E(x,y) = {xy\W\xy) = ^W ijk ix*y*x k yi > 0, (11) 

ijkl 

where {xt, yi} are the components of |x), \y) in some ba- 
sis, and Wijki are the matrix elements of W in the same 
basis. Equation (11) states that the biquadratic hermi- 
tian form E associated with W must be positive semidef- 
inite (PSD). It is not hard to show that all of the EWs 
generated by Eqn. (10) satisfy the relation 

{xyx\Z E w ® I\xyx) = (xyx\(Z + Z x + Z 2 )\xyx) 

= (xyx\Zo\xyx) + (x*yx\Z'[ A \x*yx) 
+ {xy*x\Z^ B \xy*x). (12) 

Sincc Zf),Z\ A and Zj B are positive by construction the 
biquadratic hermitian form E(x,y)(x\x) has a decom- 
position as a sum of squared magnitudes (SOS). This 
guarantees that E(x,y) is PSD. It can be shown that 
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our first separability criterion detects all entangled states 
that possess an EW such that E may be written in this 
forní. The dual program to our initial SDP may be 
interpreted as a search for an entanglcment witness of 
this kind. Equally, the Peres-Horodecki criterion detects 
thc entanglcment of those states that possess entanglc- 
ment witnesses for which (11) may be written directly as 
a SOS — the decomposable entanglement witnesses [16] 
such that W = P + Q Ta for some PSD P and Q. In gen- 
eral, if there is no EW W such that (11) is a SOS, we can 
search over W for which (11) is a SOS when multiplied 
by (a;|a;} fe_1 for some k,l > 1. By duality, this 

corresponds to our (k, l) separability criterion. 

As an example illustrating the methodology, consider 
the state described in [3, Section 4.6], given by: 

2 , , , , ce 5 — a „ „ , 

Pa = ^+)W+\ + ^+ + ^j-P<7+P, (13) 

with < ce < 5, |^ + ) = ^ Elo l«>. ct + = |(|01>(01| + 
|12)(12| + |20)(20|). Notice that p a is invariant under 
the simultaneous change of ce — > 5 — ce and interchange 
of the parties. The state is separable for 2 < ce < 3 and 
not PPT for ce > 4 and ce < 1 . Numcrically entanglement 
witnesses could be constructed for p a in the range 3 + e < 
a < 4 (and 1 < ce < 2 - e) with e > 1CT 8 . A witness for 
ce > 3 can be extracted from these by inspection: 

Z EW = 2(|00)(00| + |11)(11| + |22)(22|) + 

+ |02)(02| + |10)(10| + |21)(21| - 3IV+XV+I- 

This observable is nonnegative on separable states: 

2(xy\Z EW \xy)(x\x) = \2x xiy^ - x 2 x y{ - XiX 2 y^\ 2 
+ \2x x* y - 2x 1 xf ) y 1 + xixly - x 2 x* y 2 \ 2 
+ \2x xly 2 - 2x!X 2 yi + x 2 x* 2 y 2 - x x* 2 y \ 2 
+ \2x x\y - 2x 2 x* 2 yi + x 2 x\y 2 - XYx\y^ 
+3 \x 2 x yl - xix 2 y^\ 2 + 3 \xix*y - x 2 x* y 2 \ 2 
+3 \x 2 x* 2 y 2 - x x* 2 y {) \ 2 + 3 \x 2 x* iy2 - x^ly^ 2 > 0. 

The expected value on the original state is Tr[ZEwPa] — 
j(3 — ce), demonstrating entanglcment for all a > 3. 

The reformulation of our separability tests as a search 
for SOS decompositions of thc forms E(x,y) provides 
connections with existing results in real àlgebra (see [17] 
for a discussion of the SDP-based approach in a general 
setting). By Artin's positive solution to Hilbert's 17th 
problem, for any real PSD form /(x) there exists a SOS 
form h(x), such that the product /(x)/i(x) is SOS [18]. 
Finding such an h(x) and SOS decomposition proves that 
/ is PSD. For a fixed SOS form h(x,y), we may write a 
SDP that attempts to find EWs such that h(x, y)E(x, y) 
is SOS. In our hierarchy of criteria the form h is restricted 
to be (x\x) k ~ 1 {yly)^ 1 . While it is conceivable that every 
PSD bihermitian form is SOS when multiplied by appro- 
priate factors of this kind, currently we do not have a 
proof. It is known that deciding whether a form is posi- 
tive is NP-hard and so this connection to positive forms 



also promises to shed light on the computational com- 
plcxity of the separability problem. 

In this letter we introduced a hierarchy of separabil- 
ity tests that are computationally tractable and strictly 
stronger than the PPT criterion. Only the second step in 
this sequence of tests was required to detect the entan- 
glement of a wide class of known PPT entangled states. 
The method is based on the application of semidefinite 
programs. By exploiting the duality property of these 
problems, we showed how to construct entanglcment wit- 
nesses for states that fail any separability test in the se- 
quence. These numerical results can also be very help- 
ful in finding analytical expressions for the entanglement 
witness. The application of this approach to the charac- 
terization of positive maps will be reported elsewhere. Fi- 
nally, the wide range of applications of semidefinite pro- 
gramming, along with the work reported here and in [19], 
suggests that it may become a useful tool in quantum in- 
formation and in quantum theory in general. 
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